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Abstract 

The Grover walk, which is related to the Grover's search algorithm on a quantum computer, is 
one of the typical discrete time quantum walks. However, a localization of the two-dimensional 
Grover walk starting from a fixed point is striking different from other types of quantum walks. The 
present paper explains the reason why the walker who moves according to the degree-four Grover's 
operator can remain at the starting point with a high probability. It is shown that the key factor 
for the localization is due to the degeneration of eigenvalues of the time evolution operator. In 
fact, the global time evolution of the quantum walk on a large lattice is mainly determined by the 
degree of degeneration. The dependence of the localization on the initial state is also considered 
by calculating the wave function analytically. 
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I. INTRODUCTION 

The quantum walks are roughly classified into discrete time quantum walks 
and continuous time quantum walks 0, Q • We focus on the discrete time quantum walks on 
a square lattice. The study of the discrete time quantum walks was begun by Aharonov et 
al. l| in the early 1990s, then it has been investigated by a number of groups. The discrete 
time quantum walk evolves by repeating simple quantum operations, and it is expected to 
be realized in a quantum computer. The Grover's search algorithm 9], which is one of the 
most famous quantum algorithms, is especially related to a discrete quantum walk lOl lll|. 
Recently Shenvi et al. [12] actually proved that a discrete, coined quantum walk can equal 
Grover's algorithm. For an introduction of the implementation by a quantum computer, see 
Travaglione and Milburn 0], for example. 

The recent concentrated studies make clear mathematical properties of the one- 
dimensional quantum walks. In particular the one-dimensional Hadamard walk is studied in 

n nn nnn 

detail [H IIS lla 1121 lla 113 ■ In contrast with one-dimensional quantum walks, little about 
high dimensional quantum walks is known Thus the purpose of this 

study is to investigate a two-dimensional quantum walk called "Grover walk" . A pionee ring 
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work for the Grover walk was done by Mackay et al. |20j. Very recently Tregenna et al. 
showed numerically that the quantum walker who is controlled by the Grover's operator is 
observed at an initial location with a high probability. In this paper, this phenomenon is 
referred to as "localization". They showed also that the quantum walker starting from a 
special initial state spreads out numerically. 

The first question we have to ask here is whether the localization remains even after a 
sufficiently large time. Unfortunately numerical simulations can not give us the exact answer 
on this problem. Hence it follows that we have to calculate the wave function rigorously. 
Secondly we ask why the localization is observed only in the Grover walk. There are many 
different quantum walks, however, the localization is not observed except the Grover walk 
in our knowledge. Thirdly we would like to know the dependence of the localization on the 
initial state. We will answer these questions in following sections. 

The paper is organized as follows. After defining the Grover walk in section H, we calculate 
eigenvalues and eigenvectors of the time evolution operator to obtain the wave function in 
section HI. Section IV treats the wave function at the origin and the time-averaged probability. 
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Using the results we show that localization remains even if the system size is infinity. In 
section V, we concentrate our attention to the localization on an infinite lattice and explain 
the reason why the Grover walk is special. Furthermore we consider the dependence of the 
localization on the initial sate and show that the localization disappears at a certain initial 
state. 



II. DEFINITION OF THE TWO-DIMENSIONAL QUANTUM WALKS 
A. Time evolution of the two-dimensional quantum walks 

The Grover walk considered here is a kind of discrete time quantum walks. Thus we begin 
with defining the two-dimensional quantum walk on the square lattice Z N = {(x, y) G Z 2 \ — 
(N-l)/2 <x< (JV-l)/2, -(JV-l)/2 < y < (N-l)/2} with periodic boundary condition. 
In this paper we assume that the system size N is odd. There are four quantum states at 
each site: "R" , "L" , "U" and "D" corresponding to right, left, up and down, respectively. The 
value of wave function for one of these states S G {R, L,U, D} at the position (x,y) and 
time t is written by \S,x,y,t). The time evolution of \S,x,y,t) is determined as follows: 

\R,x,y,t + 1) = a u \R,x - l,y,t) + a 12 \L,x - l,y,t) + a 13 \U,x- l,y,t) + a u \D,x- l,y,t), 
\L,x,y,t+ 1) = a 2 i\R,x + l,y,t) + a 22 \L,x + l,y,t) + a 23 \U,x + l,y,t) + a 24 \D,x + l,y,t), 
\U,x,y,t + 1) = a 31 \R,x,y- l,t) + a 32 \L,x,y- l,t) + a 33 \U,x,y- l,t) + a u \D,x,y- l,t), 
\D,x,y,t+ 1) = a,4i\R,x,y + l,t) + a 42 |L, x,y + l,t) + a 43 \U,x,y + l,t) + au\D,x,y + M).(l) 

This evolution is characterized by the next matrix: 



A = 



a ll °12 °13 a 14 

a 2 \ a 22 a 23 a 2 ^ 

031 °32 O33 034 
°41 a 42 °43 °44 

The matrix corresponding to the Grover walk is defined by 



(2) 



A. 



(3) 
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We introduce here other two-dimensional quantum walks to compare with the Grover walk 
by setting following matrices: 



A, 



A 

u 


U 


1 

V2 


1 

V2 


n 
U 


n 
U 


1 

V2 


1 

V2 


i 

V~2 


l 

V2 


I) 


(j 


1 

v/2 


1 

V2 





o_ 


1 

V3 





1 


i 

V3 





1 


1 

V3 


1 

V3 


1 


1 

V3 


1 

V3 





1 

V3 


1 





1 

V3 J 



(4) 



(5) 



We define the wave function of the total state at time t by 
1>(t) = (V'o,o(*),V'i,o(*) j ^,o(*) j ---,V'jv-i,o(*), 

where T means the transposed operator and 

^x, y (t) = (\R,x,y,t), \L,x } y,t), \U,x,y,t} } \D,x,y,t)) T . 



(6) 
(7) 



(8) 



If the initial state ip(0) is given, then the wave function i/j(t) is calculated by the iteration 
The iteration can be expressed more compactly by introducing a 4N x AN unitary 
matrix M satisfying ip(t + 1) = Mip(t). 

The probability of observing the quantum walker at a given point (x, y) and time t starting 
from an initial state ip(0) is defined by 



P(x, y, t; ip(0)) = ^ ( s > x ' V' t\ S , x , Vi 

Se{R,L,U,D} 



(9) 



Fig. 1 shows the probability distribution P = P(x, y, t; ^(0)) corresponding to matrices (a) 
A (Grover walk), (b) Ai and (c) A 2 at t = 31 on the lattice with N = 51 starting from a 
pure initial state ip(0) = (1,0,..., 0) T . In contrast with (b) and (c) localization at 

the origin can be seen only the Grover walk 
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III. EIGNENVALUES AND EIGENVECTORS OF THE MATRIX M 



A. Eigenvalues 



To express the wave function as a function of t explicitly we consider the eigenvalues and 
eigenvectors of the matrix M. These are easily obtained by the using the Fourier transform. 
According to the previous studies jq, the eigenvaleus of the matrix M are given by a set of 
eigenvalues of the following matrix 



{A) 



u~ n 

u n 

uj- m 

uj rn 



0-11 a>Y2 Oi3 Cti4 

021 C-22 «23 &24 

«31 a 32 «33 O34 

O4I Cf-42 O43 (I44 



(10) 



where to = e 2m / N , The integers n and m are quantum numbers in a wave number space and 
they take values between and N — 1. Since there are four components in H nrn ^(A), the 
number of eigenvalues is (4iV) 2 , if not consider the degeneration of eigenvalues, and each 
eigenvalue is labeled by n,m e {0, 1, . . . , N — 1} and k e {1,2,3,4}. As a result, when 
n m, the eigenvalues of the matrix M corresponding to the Grover walker, \n,m,ki are 
given by 



A n ,m,2 = 1) 



- cos £ m - cos £ n - \ -4 + (cos £ m + cos £ n ) 



m,3 



A 



- cos £ m - cos i n + \ -4 + (cos £ m + COS f „ 



n,m,4 



(11) 

(12) 
(13) 

(14) 



where £j = 2jn/N. When n — m, the eigenvalues are written as 



An,n,l 
A n ,n,2 
An, n, 3 
An,n,4 



-1, 
1, 

— U)~ 



(15) 
(16) 
(17) 
(18) 
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B. Eigenvectors 



We write the eigenvectors corresponding to the eigenvalues \ n ,m,k as 4> n ,m,k and we let 
4>i,n,m,k be i-th element of n , m , fc . If -(iV - l)/2 < x, y < (N - l)/2 and 1 < j < 4, we find 
integers x, y and j satisfying an equation ANy + Ax + j + 2N 2 — 2 = i for a given natural 
number i. Using these x, y and j, the i-th element of (f) n ,m,k, which generates an orthonormal 
basis is given by 

.nx+my 



H,n,m,k 



k^ 



N \ J2j=l \ V j,n,m,k\ 



(19) 



where t>j, n ,m,fc is the j-th element of the eigenvector of H n ^ k . We show a set of eigenvectors 
of v n ^ k = (ui,„, m ,fc, • • • , w 4 ,„, m ,fe) T in the following: 
Case 1: m = 0, n > 0, k = 1 



v n ,o,i 



1 

-1 





Case 2: m = 0, n > and k > 1 



^n,0,fc + Ai 
^n,0,fc + /3n 

2 ^n,0,fc + PnK,0,k - Ai _ 



Case 3: n = m 



(20) 



(21) 



K,n,fc|l < < 4} = < 
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(22) 



Case 4: n + m = N 



6 



{v^N-nA 1 <k<4}=< 
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4 



(23) 



Case 5: otherwise 



^n,m,k 



^n,m,k + + Pm)^n,m,k + Oi n (3 m 

CX-nfim^n m k {Pm ^rt'^'") An.m.fc "I - Confirm, 

+ (1 + a£ + a n (3 m )\ 2 n m ^ - a n (3 n 



(24) 



where av 



a; 



and /3„ 



a; 



IV. WAVE FUNCTION OF THE GROVER WALK 

A. Expansion of wave function in terms of eigenvalues 

Since we have obtained the complete eigenvalues and normalized orthogonal eigenvectors, 
we can express the wave function as a function of time. Before we present the wave function, 
we number the states R, L, U, and D from 1 to 4, respectively. Let l(S) be the number of 
the state "S". Then we have the wave function \S,x,y,t): 

4N 2 N-lN-l 4 

\s,x,y,t) = y^</> 

i,n,m,h ( f ) j t n,m,k 

j=l n=0 m=0 k=l 

where % = 4Ny + 4x + l(S) + 2N 2 — 2 and if)j(0) is the j-th element of the initial vector ip(0). 

As shown in (|14p. the eigenvalues degenerate strongly. Thus we try to expand the wave 
function by distinct eigenvalues. The eigenvalue X n ,m,i = — 1 always exists for any combina- 
tion n and m. Furthermore the A nm 3 and A nmj 4 become —1 for n = m = 0. If n > and 
m > 0, then the eigenvalues \ n , m ,k are distinct for fixed n and m. Therefore the condition 
^n,m,k = ^n',m',k is equivalent to the following condition 

cos £ m + cos £ n = cos £ m / + cos . (26) 
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A set of a pair (n 1 , m!) belonging to the same eigenvalue X n , m ,k for /c > 2 and n > is given 
by 



Q(n, m) = < 



{(n,0),(0,n),(N -n,0),(0,N -n)} for m = 0, 
{(n,n), (n, N — n)} for n — m, 
{{n, to), (n, N — m) , (N — n, m) , (N — n, N — m) , 
(to, n), (to, N — n), (N — m,n), (N — to, N — n)} otherwise. 



To express the wave function compactly using the set f2(n, m) we define the following func- 
tions: 



c i,j,0,0,k 



H,j,n,m,k 



Vi,o,o,kVjfl,o,M Q ) 



\JYa=1 \ V i,0fl,k\ 2 \jYl)=l Ko,0,fc| 2 

V i,n',m',kV* tn ^ ml ^M 

n',m'£n(n,m) \/J2t=l \ V i,n',m',k \ 2 \ Y^j=l \ v j',n',m',k 



(27) 



for n > and to > 0. 



(28) 



We note here the reason why the system size is restricted to odd in this paper. In the case of 
odd, the degree of degeneration of eigenvalues is eight at the most except for the eigenvalues 
— 1 and 1. On the other, the eigenvalues with large degree of degeneration exist in addition 
to the 1 and —1 in the case of even. This fact does not lead us to fatal difficulty, but the 
calculation becomes more complicated than that in the odd case. 
We now have another expression of wave function: 



4JV 2 

\S,x,y,t) = 



N-l 
2 4 



Ci,j,i + Cij-i(-iy + ^2 ^2 c M>Afc^n,o,fc 

n=l k=3 

», , . V-3 V-l 

JV-l 4 ~2~ ~2~ 4 

+ Ci >h n > n > k ^n,n,k + ^2 Ci d,n,m,k^n jm ,k 



n=l k=3 



n=l m=n+l k=3 



(29) 



where 



N-l 
2 



JV-l 



JV-3 N-l 
2 2 



C it j,l - Ci,j,0,0,2 + ^2 C M>A2 + Ci '0,n,n,2 + ^ ^ Qj,n,m,2, 



(30) 



n=l 



n=l 



JV-l 
2 



71=1 m=n+l 
JV-l 



N-3 iV-1 
2 2 



- Qj',0,0,1 + Qj'0,0,3 + Q,j,0,0,4 + ^ c ij,n,0,l + ^ C M>,n,l + C M>,™-l( 31 ) 



n=l 



n=l 



n=l m=n+l 
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and i = ANy + Ax + l(S) + 2N 2 - 2. In the formula (JHJ), C i>jtl and Cy^-l)* give the 
contribution of the eigenvalue 1 and —1 to the wave function. The remaining terms are 
corresponding to the eigenvalues for k = 3 and 4. 

We move on more specific cases. In order to show the localization we calculate \R,0,0,t) 
for the Grover walk starting from a pure state \R, 0,0,0) = 1. We refer this special initial 
state as 4>r. The values in (|2~7jl . (|2*5j) and J3UJ), (JHH can be obtained after some algebra by 

N 2 1 iV 2 

= — , Ci,i -l = 2 + — ' ( 32 ) 

Cl,l,n,0,fe = 1) for v ^ ( 33 ) 
Cl,l, n ,n,l — c l,l,n,n,2 = C l,l,n,n,3 = 0, C^o.n/M = 1, (34) 
Cl,l, n ,m,fc = 2, for Vfc. (35) 



V. TIME- AVERAGED PROBABILITY 

A. Definition of time-averaged probability 

We begin with considering the probability P(S, t; O , ^V") that a walker in the state "S" at 
the origin on a square lattice with size N starting form an initial state 0o- The probability 
P(S, t; 0o, N) is calculated from the relation 

P(S, t; 0o, iV) = (5, 0, 0, t\S, 0, 0, t). (36) 

The probability P(5*, t; 0o> N) dose not converge to a fixed value in the limit t —>■ oo in con- 
trast with classical random walks. Thus we introduce time-averaged probability P(S; 0o, N) 
defined by 

1 

P(S; 0o, iV) = lim - P ( S > *; ^o, N). (37) 

t=0 

Let us calculate P(S; 4>r, N), which is the time-averaged probability starting a pure state 
<Pr. Submitting the wave function with coefficients (|32 |1 -(|35 )1 into the definition (jHEJ), we find 
the cross terms in the form \ n ,m,k^n> ,m' ,k> ■ If the eigenvalue X n ,m,k is different from X n ' >m ' t k', 
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the time-averaged value lim^oo J2t=o (^n,m,kY m , \Sf jT vanishes. Thus we have 

I ■' V : ~ I \ — ■ : 1 I : \ : _ ~" 

m 00, n) = i 



4AT 2 
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4Af2 


2 v-i 


4AT2 
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2^^,3,-1 


+ EE 


^ ^ Ci,j,nfi,k 


. J ' =1 




3=1 


n=l fe=3 





JV-l 4 

+EE 

n=l fc=3 



4Af2 

^ Ci,j,n,n,h 
3=1 



N-3 JV-l 
2 2 



+ E E E 

n=l m=n+l fc=3 



AN 2 

^ ^ Ci,j,n,m,k 
3=1 



(38) 



Plugging equations ([3*2 |) -([3l) ]) into (jSHJ), the time-averaged probability P(R; (pR, N) is given 
by 

5 



Pfifc R , AO = - + — - A 



4AT4' 



(39) 



The time-averaged probability ^(S 1 ; N) is a monotone decreasing function in N and it 
converges to 1/8 in the limit N — > oo. Thus P(S; 4>r, N) is larger than 1/8 for any odd N. 
It means that the quantum walker centralizes at the origin. 

We must pay attention to the dependence of the wave function on the parity of time. 
The value of wave function at odd time is small in comparison with the value at even time. 
It is similar to the fact that the probability of return to the origin at odd time is zero in a 
classical random walk on an infinite square lattice. Let P e {S; 4>r, N) and P (S; <J)r, N) be 
the time-averaged probabilities over even time and odd time, respectively. Then we have 

13 2 5 



Pe(R](f>R,N) 

P (R;cf> R ,N) 



4 2N 2 N 3 4iV 4 ' 
1 2 5 

Jp ~~ iV3 + AN 1 ' 



(40) 
(41) 



The probability averaged over odd time P (R; 4>r, N) converges zero in the limit iV — > oo. 
Thus there is a relation P e (R; 4>r, N) = 2P(R; <J)r, N) in the limit iV — > oo. 



B. The time-averaged probability in the limit N —* oo 

We now proceed to the probability P(S,t; 4>r, N) in the limit iV — > oo. The only first 
and second terms in (JHHJ) remain in the limit of N — > oo. Thus its calculation becomes easy 
as shown below: 

PooOS-jfo) = lim P(S;<p R ,N), (42) 
= ^j-A^if + lCi^l 2 ). (43) 
10 



Furthermore the constant and the single summation in and Cij,-i do not contribute 

to Poo(S; 4>r)- Consequently, the probability P^ (S; 4>r) is given by 



JV-3 V-l \ 1 / N-3 V-l 

2 2 \ / i 2 2 



P^S; ( j )R )= llun^Yl E c ^.--i + j^^E E c ^n, m , 2 • (44) 

V n=l m=n+l / y n=l m=n+l / 

Let consider time-averaged probabilities in N — > oo for all possible states at the origin. 
The value of Cj j0 ,n,m,i and Cj i0 ,n,m,2 as function n and m are given by 

2 (cos£ + cos£ 2 cos£ 

in 

COs£ n ) 

C2,0,n,m,l — Z~~ 7 \ 7 1 V^S 

-2 + cos 4 m + cos £ n 
2 (cos£ m + cos£ n + 2 cos£ m cos£ n ) 

C2,0,n,m,2 — — Z j Z , (40 ) 

2 + COS 4 m + COS £ n 

8sin 2 (e m /2) sin 2 (q/2) . 

Ci,0,n,m,l = ~ZZ Z 7 for 1 = 3 ' 4 ' 4T ) 

2 - cos 4 m - cos £ n 

8cos 2 (e m /2)cos 2 (q/2) 
Cj,o,n,m,2 = ~ : t — ; t for 2 = 3,4. (48) 

2 + COS 4m + COS £ n 

The double summations in f)44j) in the limit iV — > oo are calculated by replacing the sum- 
mations into the following integrals: 

y-3 v-i 



i^iv^E E c2 a^,i 



iV^oo A/ z -t— ' ' o-k 1 \ n In — 2 + COS x + cos -u 

n=l m=n+l 



2 2 

i^ooiV^E E C 3An, m ,l 
n=l m=n+l 



1 






/ ay- 

'0 


8tt 2 


JO 


1 


1 






4 ~ 


1 

7T 






1 

27T 2 


7T 

/' 

Jo 


fix 


/ 


1 


1 






4 ~ 


27r" 







(49) 



JV-3 jV-1 

8 sin x sin y 



(50) 

Similarly we can calculate the summation for C2,o,n,m,2 and C3 t o,n,m,2j and we have the same 
values 1/4 — 1/vt and 1/4 — 1/27T. Finally the time-averaged probabilities are given by 

Poo(R-Ar) = I, (51) 
Poo{L-Ar) = 1 + ^2--, (52) 

8 7T Z 7T 

Poo{U\<t> R ) = - + At-— , (53) 
111 

^oo(-D; <Pr) = - + — ? • (54) 
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Summing P 00 (S'; 4>r) over all possible states, the time-averaged probability that a quantum 
walker exists at the origin is 



PMr) = Poo(R; 4>R) + Pvo{L\ <Pr) +Pao(U; (j) R ) + Pvo{D\ 4>R 

1 3 2 

2 f 2 f 



o+^-T (55) 



C. Dependence of the time- averaged probability on an initial state 

If the Grover walk starts from the pure state, then the localization of the quantum walker 
at the origin was shown in the previous section. The numerical results obtained by Tregenna 
et al. (3], however, suggests that the time-averaged probability at origin of the Grover walk 
with a certain mixed initial state becomes zero. To confirm this observation we consider the 
wave function at the origin starting the next mixed initial state 

<f>(a, (3) = a\R, 0, 0, 0) + /3\L, 0, 0, 0), (56) 

where |a| 2 + |/3| 2 = l. In the formula in (}2*9"|) the coefficients corresponding to j > 2 are zero 
and the coefficients corresponding to j = 1 are already computed in (|35|). Thus we consider 
the only coefficients corresponding to j = 2. Suppose that ^(0) = 4>l = (0, 1, 0, . . . , 0) T . 
Then we obtain 

N-1 

c - 7 + 3N V " 

^1,2,1 - -J + -7T ~ 2^ 



^ 3 + cos £ n 



N-3 N-1 

2 (cos £ m + cos £ n + 2 cos £ m cos £ n ) 



hrk+, 2 + 008^ + 008^ 



N-3 JV-l 
2 2 



3 ^ , ^ 2 (cos£ m + cos£ n - 2 cos£ m cos£ n ) 
Cl ' 2 '- 1 " 4 ~ T + ^ \ -2 + cos^ + cos^ ' (58) 

n=l m=n+l 

4 

c l,2,n,0,3 — c l,2,n,0,4 — ~ 1 + — — , (59) 

3 + COS i n 

Cl,2,n,n,3 = Cl,2,n,n,4 = 0, (60) 

4(cOs£ m - COs£ n ) 2 

Cl,2,n,m,3 = Ci i2 ,n,m,4 = 7 77^ 7^F~\ a C T ^ 111 > H ' ( 61 ) 

6 - cos(2£ TO ) - cos(2£ n ) - 4 cos4„ cos £ m 

We now can calculate the time-averaged probability P(R; 4>(a, (3), N) for any odd system 
size N by combining (J32J) - f)35[) with ()57|) - (|61|) . However it is rather complicated. Thus we 
show only the result in the limit of iV — > oo. The time-averaged probability P^R; <fi(a,j3)) 

12 



is obtained by replacing the summations to the integrals in the same way described in the 
previous section, and it becomes 



P^R-^a,?)) 



Similarly we have the time-averaged probability corresponding to the "L" state 

2 



Poo{L; (f){a,pj) = - 



(3+[l--)a 



71 



(62) 



(63) 



Fig. 2 shows Pr((x) = Poo(R', 4>{oc, yl — a 2 )) and Pl(oc) = Poo(L; <p(a, y/1 — a 2 )) for a G 
[—1, 1]. The time-averaged probability Pr{o) becomes zero at a m i n given by 



a, 



1 - 



7T 



16-87T + 27T 2 ' 

and it takes the maximum value at a max given by 

7T 



(64) 



(65) 



V16 - 8tt + 2vr 2 

Since the value a = 1 and a = — 1 indicate the pure state, it is found that the time-averaged 
probability Pp( a) takes the maximum at the mixed initial state. As shown in the numerical 
calculation the quantum walk starting from a certain initial state spreads out. Al- 
though the time-averaged probability Pr{o) becomes zero at a = a m i n , the quantum walker 
remains at the origin. Put another way, the component of the time-averaged probability 
corresponding to the "R" state converges to zero, but other component remains positive. 
Therefore we can say that the quantum walker in the state "R" perfectly converts into other 
state in the limit t — ► oo at a = a m i n . 

Let us consider the time-averaged probability for more general initial states defined by 

<j>(a, (3, 7 , C) = a\R, 0, 0, 0) + (3\L, 0, 0, 0) + 7 \U, 0, 0, 0) + (\D, 0, 0, 0), (66) 

where | or | 2 + |/5| 2 + |7| 2 + \(\ 2 = 1. Repeating the procedure described above yields the 
time-averaged probability for this initial state: 

a 



1 2 l a 1 



S 2tt 2 2tt^ 7 + ^ 



2\/2 V 8 7T" 7T 

Taking symmetry into account, the other components are given by 

P^L-Aia,^ 7,0) = Poo(i?;0(A«,7,O), 
P^U-Aia, (3,^,0) = Poo(R; 0(7, C, «,/?)), 
Pao(D; <p(a, (3, 7, ()) = P^R; 0(C, 7, «, 13)). 



(67) 

(68) 
(69) 
(70) 
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The condition for which all components in (|bT|) - (|70p become zero is obtained by 

a = e ie /2, (3 = a, 7 = -a, ( = -a. (71) 

Setting 9 = 0, we confirm the numerical results given by Tregenna et al. [22], which 
claim that the probability of Poo(S; <fi(a, j3, 7, ()) becomes zero for any state 5* by setting 
a = (3 = —7 = — ( = 1/2. Fig. 3 shows a snapshot after 30 steps on a square lattice with 
N = 51 starting from an initial state with 9 = 1/3 in (|71j) . One definitely finds that the 
localization disappears. On the other hand, the summation of Poo(S; j3, 7, ()) over all 
possible states takes a maximum values 2 + 8/tt 2 — 8/7r=0. 26409.. ata = /? = 7 = £ = 1/2. 
Thus the time-averaged probability over even time is larger than 1/2. 

D. Conditions of the localization 

We saw that a spike exists at the origin in Grover walk, but the two-dimensional quantum 
walkers governed by the matrix Ai and A 2 spread out. The significant difference between 
the Grover walk and other quantum walks is the degree of eigenvalues. For example, the 
eigenvalues of H n rn with A x are given by 



A n 


m,l 




a/ i COS £ n cos £ m 


"I - fn,mi 


A n 


m,2 




- y/i COS £„ COS Cm 


fn,mi 


An 


m,3 




y/i cos £ n COS £ m 


fn,mi 


_ A n 


m,4 




-■s/i cos £ n cos £ m 


fn,m} 



(72) 



where 

/n,m = \J sin 2 £ m + cos 2 f n cos 2 £ m . (73) 

We find no common eigenvalues to all value n and m such as —1 and 1 in \ n . m ,k- If a quantum 
walker exists only at the origin initially, the coefficient Cjj^m^ in (|38j) takes non-zero value 
for only one j. Therefore the summation over j in (|3*Hj) does not depend on the system size 
and the summation over n and k increases with the system size. Accordingly, the order of 
third and fourth term in (|38|) is 0(N~ 3 ). Similarly the order of the fifth term in (J38)) is 
0(N~ 2 ). As a result, these terms vanish in the limit iV — > 00. As shown in ()30j) and (|3*T| . 
since C^i and Cij-i contain double summations, the order of them is 0(N 2 ). This large 
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contribution to the time- averaged probability comes from the fact that the eigenvalues — 1 
and 1 exist for any n and m in common. From the above consideration one conjectures 
that there is another two-dimensional quantum walk in which the walker centralize at the 
origin. Thus we find a new matrix A 4 such that if nj , m (A 4 ) contains eigenvalue —1 and 1 
independently on the value n and m. Suppose that the matrix A 4 is real symmetric matrix. 
Then a possible matrix is given by 



Aa 



-p q y/pq 

q —p y/pq y/pq 

pq y/pq —q p 

pq Jpq p —q 



(74) 



where q — 1 — p. The Grover walk is corresponding to p — 1/2. As an example we show the 
probability P(x,y,t; <Pr) of a quantum walk which is characterized by the matrix A 4 with 
p = 1/3 and q = 2/3 in Fig.4. We clearly recognize the existence of a spike at the origin. 



VI. SUMMARY 



We have shown analytically that the localization of the Grover walk at the initial position 
can be surely measured for any odd system size. On the other hand, we have also shown that 
we can chose special initial states with which the quantum walk disappears at the initial 
position in N — > 00. As pointed out by Tregenna et al. this different behavior can 

be used to control the Grover's search. We here summarize the reason why the localization 
exists in the Grover walk starting from a local position. 

There are 16A 2 quantum states in the Grover walk on the square lattice including N 2 
sites. Therefore the number of eigenvalues and eigenvectors of the time-evolution operator 
is also 16 A 2 . Since the wave function at time t is expressed by a liner combination of the 
t-th power of the eigenvalue in which the coefficients are given by the product of component 
of eigenvectors. Each component of the eigenvector decrease in inverse ratio to the system 
size and the coefficients decrease in proportion to N~ 2 . Thus, if the quantum walk exists 
initially at an fixed point and all eigenvalues are distinct, then the probability of observing 
the quantum walk at the fixed point goes to zero by taking the system size infinity. For 
this reason, the degeneration of eigenvalues is necessary for the localization. In the case of 
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quantum walks on a circle including odd sites, the eigenvalues are distinct except for trivial 
cases Thus the localization is not observed. 

The eigenvalues of the Grover walk include —1 and 1, and the degree of the degeneration 
of them are N 2 + 2 and N 2 , respectively. As mentioned above, the each coefficient itself 
decrease in the form N~ 2 , but the degree of the degeneration is proportion N 2 . As a 
consequence, the eigenvalue —1 and 1 can positively contribute to the probability even if 
N — > oo except for the case that both coefficients of —1 and 1 become zero. If we choose 
such initial states as both coefficients of —1 and 1 corresponding to the wave function at the 
origin are zero for all states, then the quantum walker spreads from the origin. 
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FIG. 1: Snapshots of probability distribution P(x, y, t; ^(O)) at time t = 30 on a sure lattice with 
system size N = 51 with ^(0) = (1,0,..., 0) T . Each of evolutions are determined by the matrices 
Aq (Grover walk), Ay and A2 defined in ©-© corresponding to (a), (b) and (c). The central peak 
in (a) shows the localization of the Grover walk. 
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FIG. 3: The probability P(x,y,t; ^(0)) of Grover walk at t = 30 starting a mixed state ^(O) 
(e 1 / 3 , e*/ 3 , — e i//3 , — e i//3 , 0, . . . , 0) T with JV = 51. The central localization in Fig. 1 (a) disappears. 
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FIG. 4: The probability P(x, y, t; 4>r) of the quantum walk whose wave function is determined by 
the matrix A4 with p = 1/3 and q = 2/3. A sharp distribution at the origin is observed similarly 
as in the case of the Grover walk. 
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